Prom Hubbard model to t-J-U model: a canonical transformation formalism, the 
metal-insulator transition and mean-field state 
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We prove that the t-J-U model can be deduced from the Hubbard model at a large but finite U 
by a canonical transformation. We argue that the system may have a metal-insulator transition 
at a critical on-site Coulomb interaction whose value, however, is smaller than that in previous 
calculations in which the kinetic energy has a double counting. In a mean field theory and a 
special choice of the parameters, we show that the metallic state may be equivalent to the gossamer 
superconducting state proposed by Laughlin recently. 



PACS numbers:74.20.-z,74.20.Mn,72.-h,71.10.Fd 

The Hubbard model is the generic model for inter- 
acting electrons in the narrow-band and strongly cor- 
related systems [1]. Especially, since the high temper- 
ature superconductivity was discovered in the cuprates, 
the Hubbard model on two-dimensional lattices as well 
as its strong coupling limit model [2] , the t- J model, have 
been extensively studied in order to understand the vari- 
ous anomalous properties of the cuprate superconductor. 
The up-date investigations, however, can not supply a 
definitive evidence to show the stable (i-wave supercon- 
ducting ground state in these strongly correlated models 
[3]. 

It was known that while the Hubbard model is weaker 
in describing the exchange correlation, the t-J model is 
weaker in the long range charge correlation. A better 
phenomenological model to include stronger correlations 
is the t-J-U model [4] . The existence of both J- and U- 
terms is very important in a possible mechanism of the 
high T c superconductivity, gossamer superconductivity, 
proposed by Laughlin recently [3,5,6]. Although includ- 
ing both J- and U- terms in the model Hamiltonian is 
physically reasonable for both the exchange interaction 
and the on-site Coulomb interaction are always relevant 
to the real systems, there is no a microscopic deduction 
of the t-J-U model from a generic model. 

Comparing with the Hubbard model, there seems to 
be a double counting in the t-J-U model to the inter- 
acting or kinetic energies due to taking both the J- and 
U- terms into consideration simultaneously. All previous 
theories on the t-J-U model did not deal with this double 
counting since the lack of a microscopic understanding. 
In this paper, we would like to show that the t-J-U model 
in fact can be derived from the origin single-band Hub- 
bard model in a large but finite U. The double counting 
in the previous phenomenological t-J-U model was in the 
kinetic energy as we will show. An easy understanding in 
fact can be seen by checking the perturbative deduction 
from the Hubbard model to the t-J model: The J-term 
is actually from the second order perturbation in t/U 
by taking the off-diagonal term for the double-occupied 
number in the kinetic term. Namely, if the t-term is de- 
composed into Tdiag + Toff-diag, the perturbation process 
transfers T ff_diag to J-term, a virtual hopping process, 



and only Tdi ag severs as the real hopping. It can also 
clearly be seen from the canonical transformation deduc- 
tion of the t-J model [7]. 

In this paper, we would like to deal with the Hub- 
bard model in a large but finite on-site Coulomb inter- 
action by using the canonical transformation. We find 
that the effective Hamiltonian can be written as the sum 
over the Hamiltonians acting on a subspace of the Hilbert 
space with a fixed double occupied number. In a partial 
Gutzwiller projection, we examine the variational ground 
of the effective Hamiltonian using the Gutzwiller approx- 
imation. A metal-insulator transition is found at a crit- 
ical U c which however is lower than that estimated by 
Zhang [5]. The overestimate of U c in the previous works 
stems from the double counting of the kinetic energy in 
the phenomenological t-J-U model. 

Our model is [/(l)-invariant due to the particle num- 
ber conservation. If we consider a mean field state with 
the symmetry breaking, we can show that the mean field 
state with a special choice of the parameters is equivalent 
to Laughlin's gossamer superconduting state [3] . 

We start from the Hubbard model where the hopping 
energy may be dependent on the occupation of sites in- 
volved [8]. Including the on-site Coulomb interaction, 
this Hubbard model reads 



H = T + V = T 



M 

i=i 



(1) 



where M is the number of the site; 



with Cj CT a spin-er electron anahilation operator at 



site i and the kinetic term is given by 
T = Th + T c i + T mix , 

Td = — ^ t tij n i& c icr C j(J n ja, 

T m ix = — ^ ] tij n i9 c ia C j cr (^ ~ n js) 



(2) 
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Except in the mean field theory, we assume = tfj = 
j.mix _ ^ f or ^ ne neares t neighbor sites and vanish other- 
wise. 

In order to define the canonical transformation, we ex- 
plain our notations. The partial Gutzwiller projection 
operator 

N/2 

n(s) = (i - g)ui) = £ $ Dp n = ( 3 ) 

i D=0 

where < g < 1 is the Gutzwiller parameters; N 
is the electron number [9], D = Y^i v i an( i Pd = 
J2{t 1 ,...a D }Wii--- L 'iD U'j( l - v i)\ is a projection operator 
which projects a state into the subspace with a fixed 
double-occupation number D. Pq = 11(0) is the full 
Gutzwiller projection operator and 11(1) = 1. For conve- 
nience, we denote Po{g) = g D PD,P m (g) = Y.D>i p o{g)- 
The first goal of this work is to construct an effective 
Halmitonian H c g and after the partial Gutzwiller projec- 
tion, the projected effective Hamiltonian is given by 

U(g)H eS U(g) = £ P D (g)H eS P D (g) 

D 

= Y,9 2D PdH cS P d , (4) 

D 

i.e., all the off-diagonal part PD'{g)Hef(PD{g) — for 
D 1 7^ D. We shall prove that the result effective Hamil- 
tonian in which all terms keep D-invariance is given by 

H cS = T h + T d + J + V, (5) 

where 

(ij) 

- 2 n i n 3\ n n + n^nnn^nji), (6) 

with Jij — J w At 2 /U for a larger U. 

The canonical transformation for the Hubbard model 
to the t-J model was introduced by Hirsh [7]. A de- 
tailed review for the canonical transformation can be 
found in Ref. [10]. Our derivation is a generalization 
of Hirsh's D = case. Notice that PoTm^Po' = 
8d> D±iPDT m i x PD±i and (5) remains D invariant, as well 
as n(g)U(g') = U(gg'), P nD {g)P nD (g') = P V M)- 
Keeping these in mind, we do a partial projection 
U(x)HH(x) with x = g 2 l N ■ For large N, x is very close 
to 1. A straightforward calculation leads to a rewriting 
of U(x)HU(x) 

U{x)HU(x) = H (x) + H^\x), 

H (x) = H diag (x) + H v D) (x), (7) 

D=2 



where 

H diag (x)= J2 P d(x)HP d (x), 

D=0 

H^(x) = P D ^(x)TP D (x) + P D (x)TP D ^{x). (8) 

The purpose of the canonical transformation is to aquire 
an effective Hamiltonian such that PqH^Pd = 

P D H$P = for D ^ 0. This H$ is defined by 

ff^e^'n^fln^e-^. (9) 

As well-known [7,10], is determined by the self- 

consistent condition iH^\x) + [Ho(x), S^] = Oandthus 
the effective Hamiltonian reads 

H^=H (x)+ l -[S^,H^(x)}. (10) 

Solving the self-consistent condition , in a large U 
is given by [7,10] 

PqH^P « P HP - lp Q HP m HP 

« P H cg P , (11) 
P ni (x)H^P Vl (x) « P m (x 2 )HP m (x 2 ). (12) 

From the middle to the right sides in (11), the three site 
processes are neglected. If the non-double occupied con- 
straint is imposed, (12) vanishes because it is related to 
the double occupation. Eq.(ll) gives rise to the common 
t-J model. However, if the double occupation is allowed, 
we have to deal with (12). In fact, one can repeats the 
canonical transformation to (12) . We would like to rc- 

quire an effective Hamiltonian H cff ' whose off-diagonal 

part PxH^Pd = P D H$P! = for D > 1. For this 
purpose, one writes 

U(x)H^n(x) = PoH^P + H (x 2 ) + H^(x 2 ), (13) 

where H (x 2 ) = P 1 (x 2 )HP 1 (x 2 ) + P m (x 2 )HP m (x 2 ). We 
do a canonical transformation and define 

^=e^ (2, n(^n W e-^, (14) 

where is required to satisfy P S^ = S&P = 

such that PoH(x)H^U(x)Po is invariant under the 
transformation and it is self-consistently determined by 
iH {2 \x 2 ) + [H (x 2 ),SW] = 0. Hence, similar to (10), 
one has 

H$ =P Q H^P + Mx 2 ) + \\S (2 \H^{x 2 )]. (15) 

(2) (2) 

Projecting H^ s ' — > H(x)H^ a 'H(x) and repeating the sim- 
ilar procedure to deduce (11) and (12), one arrives at 
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PoH^Po « P H cB P , 
P 1 (x)H^P 1 (x)^x 2 P 1 H eS P 1 , 
P m (x)H^P m (x) « P^z 3 )^^ 3 ), 



(16) 



for a large U, where the three site processes have been 
ignored. 

Repeating this procedure, we finally have 



D=0 

= U(g)H eS U(g). 



(17) 



The last equality is because H c g is D-invariant. The 
Gutzwillcr parameter is g but not x because we are doing 
the partial projection in each time canonical transforma- 
tion. Thus, we end the proof of (4) and (5). Moreover, 
we see that, in a partial Gutzwiller projection, the varia- 
tional ground state energy is given by a polynomial of the 
Guztwiller parameter g in power of 2D. The coefficient of 
g 2D -term is the ground state energy of the system with a 
fixed D. Using g as a variational parameter may be con- 
venient for the numerical simulations. Since a larger U 
will be considered, the optimal value of g is expected to 
be much less than 1. Thus, it may be enough to calculate 
only a few terms of the polynomial in the energy. The 
numerical work for this model is in progress by using a 
variational Monte Carlo simulation [11]. In this work, we 
are going to discuss an infinite system by the Gutzwillcr 
approximation [12]. The variational energy per site is 
given by 



= (y>o|n(g)H cff n(.g)|y>o) 
M(v>o|n( 5 )n( ff )|^>> 

= Ud+((T h + T d )) + (J). 



(18) 



Here d 
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is the double occupied concen- 



trate. Using d as the variational parameter, Zhang has 
discussed the metal-insulator transition at a critical U c 
for an effective Hubbard model (or the phenomenological 
t-J-U model) [5]. As we have mentioned, the Hamilto- 
nian used in his paper overestimated the kinetic energy, 
comparing with our microscopic derivation of the t-J-U 
model. Meanwhile, we will see that the critical value of 
U c estimated in [5] is higher than that we present here. 
Taking the approximation J2{ij) u i n j ~ and neglect- 
ing the viVj term in (5), the variational condition is very 
simple 



0. 



(19) 



here, J ex is the Si • Sj term in the Hamiltonian; U c s = 
U — AJn because the contribution from z/jjij term in(5) 
and g' t is modified from the original Guztwiller factor g t 
because of the absence of T m ; x . gj is the same as that in 
the projected Fermi liquid state [12]. We have, in fact, 



9 J = 



9t = 



4(1 -(5 -2d) 2 
(l + <5) 2 (l-<5) 2 ' 
2(1 -5 -2d) 

TTs 



■d d(l-J) 



1-5 (l + 5) s 



(20) 



while g t = (f^ggjnjlp , where 5 = 1 — n. At half filling, 
5 = and it is easy to see g' t = 4d(l — 2d) = gt/2. 
Hence, at half filling, the critical point of the metal- 
insulator transition is at U c = — 4(T) + 16( J ex )o + 4J. 
Using (T) « -2y/2tC and (J ex ) « -(3/4)JC 2 with 
C w 0.479 [5], U c = 5 At + 1.25J. This estimate for 
the critical point is much less than U^ h = 10. 8t — 2.75 J 
which was given by Zhang [5]. For example, if J ~ QAt, 
U c <~ 5.9t < U, while in Rcf. [5] the estimate value 
Uf 1 ~ 10.lt. If one estimates J w At 2 /U, U ~ lOt 
which is near U% but larger than U c . This means that 
for such parameters J and U, the system is in an abso- 
lute insulator state in our theory but is near the critical 
point in the phenomenological model. If J ~ 0.67i ( cor- 
responding to U ~ 6t), U c <~ 6.2t is near U, while in Ref. 
[5] the estimate value U^ h ~ 9.0t > U. In this case, the 
system is near the critical point in our theory but in an 
absolute metal state in the phenomenological model. Al- 
though the quantitative value of U c needs to be refined 
by numerical calculations, the conclusion U c < U^ h will 
not change. 

We now are going to discuss the mean field state of 
our theory. The basic idea to go this mean field state has 
been explained in our previous preprint [13], although 
there were some computational errors in that paper. Here 
we present a renewed version of the mean field state. 
Introducing two correlation functions Ay = (c^Cj-f — 

c H c jl)o,Xij = (4] c 3l + c a c Ji)o, the U(l) symmetry of 
H c s is broken by a decomposition of the four particle 
terms [13]. According to Ay and the mean field 
Hamiltonian of (5) is given by 

H M F = - ^ (*« + O 71 ^ + n 3s) + 4f n i<r n j<r) C L C j<r 
+ + J y 1} ( n ^ + + J iT n i° n jv) 

x (-inAf.c^ + Ayct-4) 

+ U n^riii - 2Jy (1 - Ajn^n^ (21) 

' to) 

- £ Jy (A|Ay | 2 + ^(1 - B)\ Xij \ 2 ){ni + n 3 ) 

m 

+ y ] n ia n jS H ^ \Xij\ n iu n ja), 

where the parameters arc given by 



J (1) - — J 



J (2) - -*J. 



t$> = -$ j -0--A)XiiJii, 



(22) 
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tff =tt j +tf j -(l-B) Xji J ij . 

A and B are the variantional parameters to be deter- 
mined. On the other hand, we write down Laughlin's 
gossamer superconducting Hamiltonian 

H G -ii R N = Y,^J**, (23) 

k 

w here \xr is reno rmalized chemical potential, E k = 
\J (ffe - ^r) 2 + A 2 , and 6 kff = n^b^ir 1 ^ for 6 kT = 
UkCkt + ^kclki and &kj. = ^kCkx — ^kcl k | annihilate the 
BCS state. Expressing explicitly (23) by the electron 
operators [13,6], we have 

H G - m N = -J2 [*« + *2 (1) (»iff + 

+ ^W^ja] + ^ Jy [1 + ^a[3(n is + n j£r ) 

- a(3n lS n 3S ]{-l) a {A\ jCia c js + A ijC ^c\ a ) 

+ U G Y,n^n il - ti G N (24) 

i 

where a = 1 — g and (3 = {1 — g)/g and 

^ (2) =E§(^-/3M)e ik - ( --^, 
fe 

^•^O—E^Vke*-^-'^, (25) 

k 

and C/ G = ^Ek^[(2/3 + /3 2 K + (2a-a 2 K], MG = 
Tf J2k E k[(2a+l)vl-u 2 k }. If we identify the t-J-U model 
to the gossamer superconducting model in the mean field 
level, one requires 

A = a(3, B = 2a(3, 

* G(1) = ~t% - (1 - afixji Jij, 

t% 2) =t* j +tt j -(l-2a0) Xj iJi j , (26) 

and [in + fi G — J(12A|A T | 2 + 8(1 - B)\ Xt \ 2 ) + n,U = 
U G + 8J(l-A). 

Although we have made a formal equivalence between 
our mean field state Hamiltonian to Laughlin gossamer 
superconducting Hamiltonian, we note that the hopping 
parameters i 6 ^ 1 ' 2 ) have run out of the practical range in 
the real materials. Thus, to show the system described by 
the t-J-U model has a gossamer superconducting phase 
described by Laughlin gossamer superconducting Hami- 
tonian, a renormalization group analysis is required. We 
do not touch this aspect in this work. However, we can 



believe there is such a superconducting phase in our the- 
ory if U < U c because the superconducting paring pa- 
rameter is determined by the optimal exchange energy 
as in the common t-J model. The renormalization of 
the hopping parameters is believed to affect the normal 
dissipation process only. 

In conclusions, we have deduced a t-J-U model from 
the Hubbard model at a large but finite U by a series of 
canonical transformations. This microscopic derivation 
showed the overestimate to the kinetic energy in previous 
work by a phenomenological t-J-U model. We then saw 
that the critical U c for the metal-insulator transition is 
much lower than the value calculated by the phenomeno- 
logical model. Through a mean field theory, we showed 
that the gossamer superconducting state may possibly 
be a symmetry breaking state of the t-J-U model as the 
BCS common superconducting state versus the Fermi liq- 
uid. Another advantage of our t-J-U model is that the 
energy was decomposed into the summation over the en- 
ergies with the fixed D. This is convenient for numerical 
simulation. 
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